Surface losses and self-pumping effects in a long 
Josephson junction — a semi-analytical approach 
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The flux-flow dynamics in a long Josephson junction is studied both analytically and numerically. 
A realistic model of the junction is considered by taking into account a nonuniform current distri- 
bution, surface losses and self-pumping effects. An approximate analytical solution of the modified 
sine-Gordon equation is derived in the form of a unidirectional dense fluxon train accompanied 
by two oppositely directed plasma waves. Next, some macroscopic time-averaged quantities are 
calculated making possible to evaluate the current-voltage characteristic of the junction. The re- 
sults obtained by the present method are compared with direct numerical simulations both for the 
current-voltage characteristics and for the loss factor modulated spatially due to the self-pumping. 
The comparison shows very good agreement for typical junction parameters but indicates also some 
limitations of the method. 
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I. INTRODUCTION 

In recent years the flux-flow (FF) dynamics in a long 
Josephson junction has attracted considerable attention 
in view of possible applications in superconducting mm- 
wave electronics^—. The FF mode appears in a junction 
immersed in a sufficiently large external magnetic field 
and can be described briefly as a unidirectional viscous 
flow of a dense train of fluxons (magnetic flux quanta). 

Recently, the FF oscillators have found various appli- 
cations, e.g. in mm-wave integrated receivers^—, however 
some important parameters, such as the radiation line- 
width are still not satisfactory. Therefore, there is still 
a need for more adequate description of a real Joseph- 
son junction operating in the FF mode, by taking into 
account some additional factors which may affect the 
current-voltage (I-V) characteristic of the junction and 
change its working conditions. 

In the majority of papers dealing with long Josephson 
junctions (see e.g. Refs. I8l4l4) a simplified version of the 
the sine-Gordon (sG) equation is considered, usually ne- 
glecting surface losses and assuming uniform distribution 
of the bias current density. Only recently, a few papers 
have been published^ - — taking into account a nonuni- 
form bias current distribution and its influence on the 
junction behavior. Moreover, in Refs. [Hill] a realistic 
model of the FF oscillator has been investigated both ex- 
perimentally and numerically, including general bound- 
ary conditions, nonuniform bias current profile, surface 
losses as well as self-pumping effects^ related to addi- 
tional tunneling of quasiparticles due to the Josephson 
radiation. 

The aim of the present paper is to present an analytical 
approach to the modified sG equation, which takes into 
account: (i) nonuniform current distribution, (ii) surface 
losses, and (iii) spatial modulation of the loss factor re- 
sulting from the self-pumping effect. Contrary to Ref. Il6l 
we assume standard open-circuit boundary conditions to 
make the influence of various effects more pronounced. 



Nevertheless, the present analysis can be easily extended 
to include also more general boundary conditions. 



In the particular case of uniform current distribution, 
the present method makes it possible to obtain fully an- 
alytical closed-form expressions describing both the su- 
perconducting phase within the junction and the I-V 
characteristic. Such a solution can be regarded as the 
first-order approximation which appears sufficiently ac- 
curate for some moderate junction parameters. How- 
ever, in the general case, particularly for very long and 
weakly damped junctions, such an analytical approxi- 
mation is only the first step in an iterative procedure, 
which has to be performed numerically. Thus, in spite 
of analytical expressions describing the superconducting 
phase within the junction, the present approach has been 
named "semi-analytical" . 



The paper is organized as follows. In Sec. II we for- 
mulate the problem, i.e. we present a generalized sG 
equation subject to open-circuit boundary conditions at 
the junction ends. Approximate analytical solutions to 
the sG equation are discussed in Sec. III. We start with 
a linearized (small-amplitude) solution and apply appro- 
priate boundary conditions. Next, some large-amplitude 
corrections are introduced, we discuss also possible self- 
pumping effects and their influence on the I-V character- 
istic. In Sec. IV analytical results are compared with di- 
rect numerical simulations. In particular, we discuss the 
influence of surface losses both for junctions of moderate 
length and for more realistic structures, such as very long 
junctions with small damping and strongly nonuniform 
current distribution. Sec. V contains concluding remarks, 
we indicate also possible extensions of the method by tak- 
ing into account more general boundary conditions. 
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II. FORMULATION OF THE PROBLEM 



III. APPROXIMATE ANALYTICAL SOLUTION 



Fluxon dynamics in a long Josephson junction is usu- 
ally described by the following modified sG equation^—: 

(f>xx - 4>tt - OL(j>t = sin <p - 7, (1) 

where <p denote the quantum phase difference across the 
barrier, a is the loss factor, and 7 is the bias current 
density. The spatial coordinate x has been normalized 
to the Josephson penetration depth Aj and the time co- 
ordinate t to the inverse plasma frequency oJq 1 , where 



A./ = (h/2fi edj l 



a/2 



CJ 



= (2ej c /HC) 1 ^ 2 , j c is the critical 



current density and C denotes the junction capacitance 
per unit area. 

However, to describe more adequately a real physical 
situation, one can consider a more general for m 16 ' 17 : 

4>xx - <fot - a(x)(f>t + f3(/)xxt = sixi(p- 7(e) , (2) 

where (3 denotes the surface loss parameter and we as- 
sume both a and 7 to be cc-dependent, taking into ac- 
count both a spatial modulation of the loss factor due 
to the self-pumping and a nonuniform distribution of the 
current density along the junction. 

For the overlap geometry one can neglect the self-fields 
and assume simple open-circuit boundary conditions 9 : 



,(±L/2) 



t (±L/2) = h, 



(3) 



where h = H ext /j c Xj, H cxt denotes the external mag- 
netic field, and L ^> 1 is the normalized junction length. 

Following Refs. ITl1 - fl4lfl8l we look for an approximate 
solution in a form of a dense fluxon train traveling on 
a rotating background. Thus, a linearized solution of 
Eq. J2j| can be written as 



<f> = (f>Q + "0, 



(4) 



where 0o = 6( x ) + is the background term (linear in 
time) and ip denotes a quasi-linear term (usually small) 
representing the motion of fiuxons and plasma waves 
within the junction. 

Substituting Eq. (gj) into Eq. Q we find 

Oxx + ^xx-^tt-aVL-a^t+P^xxt = sin(0 o + T/>)-7- (5) 

We are interested in a steady-state, time periodic so- 
lution, thus the background frequency f2 should be equal 
to the fundamental frequency of the oscillating term ip. 
Both experimental data and numerical simulations show 
that the output signal of a real FF oscillator is nearly 
sinusoidal 1,15 . Thus, it is reasonable to assume the time- 
dependence of the term ip to be harmonic in fit, and 
neglect any higher harmonics generated by the nonlin- 
ear term sin(0 o + ip)- Consequently, Eq. ([5]) can be split 
into a time-independent part and a part oscillating with 
frequency fi, while the boundary conditions © can be 
written separately for 0(x) and ip(x,t): 



x (±L/2) = h, 
iP x (±L/2)+^ xt (±L/2)=0. 



(6) 
(7) 



The method for solving Eq. ([5]) with the boundary con- 
ditions J5]) and d?]) is similar to that reported recently (see 
Ref. [H for details). Thus, below we discuss briefly only 
the main points of the analysis. 



A. Time-independent equation 

It can be easily shown that for ip oscillating with fre- 
quency Q, the nonlinear term sin^o + VO contributes also 
to the time-independent part of Eq. ([S| , and such a con- 
tribution yields, in fact, the first-order approximation of 
the I-V characteristic of the junction. Thus, the time- 
independent equation can be written as: 



a(x)fl + S(x) — 7(2:), 



(8) 



where 



S(x) = (sin(0 o + ip)) T - (1/T) / sin(0 o + iP)dt, (9) 

Jo 

and T = 2ir/n. 

Solving Eq. © for 6{x) we find 



-L/2 



[«(0n + S(0 - 7 (0]d£ 



-L/2 



dr/+Cix+C 2 . 

(10) 

Without loss of generality C2 can be set to zero while 
from the boundary condition (j6]) we find C\ — h and 



70 = a fi - 



1 



L/2 



L/2 



S(x)dx, 



(11) 



where 70, ao denote the averaged current density and 
loss factor, respectively: 



7o 



L/2 



j(x)dx, 



-L/2 



1 



L/2 



a{x)dx. (12) 



-L/2 



It follows from the relation (H|) that the time-averaged 
value of (fit, (being proportional to the constant voltage) 
is equal to the frequency i7. Thus, the expression (jTTJ) can 
be regarded as the current-voltage (I-V) characteristic 
of the Josephson junction. One can see that the total 
current density 70 consists of a linear (Ohmic) part and a 
nonlinear contribution related to the Josephson current. 



B. Time-dependent equation 

For ip sufficiently small we can write 

sin(0o + ip) — sin <j>o + ip cos 4>q. (13) 
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Noting that the term tp cos <po gives no contribution of 
frequency f2 and collecting all the time-dependent terms 
in Eq. (O we find 



Tpxx - At - a(x)ip t + /3ip xxt ~ sin0 o . 



(14) 



Since all the terms in Eq. (|T4"|) oscillate with the 
same frequency f2, it is convenient to use a complex 
notation 1 —^ 



tp(x,t) = Im ip(x)e 



jnt 



Im 



j(e+nt) 



(15) 



sm <po 

and rewrite Eq. (|14[) as an ordinary differential equation 
(1 + iftp)tp xx + (fl 2 - ia(x)fl)tp = e ie , (16) 



tp. Substituting 9 and ip into Eq. @ we obtain a new 
approximation for and consequently new approxi- 

mations for 9 and tp. 

It is clear from Eq. pT|) that 5(x) = corresponds 
simply to the Ohmic line 70 = a^Sl,, while consecutive 
iterations yield a sequence of approximations for S(x) 
describing the time-independent contribution from the 
Josephson current. 

The rate of convergence of the iterative process de- 
pends strongly on the junction parameters. As shown 
in the next Section, for moderate values of L, h, and a?o, 
the first iteration yields satisfactory results. However, for 
practical FF oscillators which are usually based on very 
long junctions, up to a few thousands of iterations are 
needed to obtain a self-consistent solution. 



iP xx + S 2 ^ = e i9 /P, 



(17) 



where P = 1 + iSl/3, S 2 = (n 2 - ia(x)Q)/P. 

The general approximate solution of Eq. (JTTJ) can be 
written in a WKB-likc forrr>2fl 



tp(x) 



1 



2iPVS 



+A e -^+B 



e -*/(») 



-if{x) 



F+{x) 



Vs 



Vs 



(18) 



where 
F ± (x) 



-L/2 



-L/2 



(19) 

The term in square brackets in Eq. (fig)) corresponds 
to a dense fluxon train moving unidirectionally along the 
junction, while the last two terms describe two plasma 
waves propagating in opposite directions. 

Using the complex formalism the boundary condition 
([7]) can be rewritten as 

M±L/2) + inp$*(±L/2) = ==> M±L/2) = 0. 

(20) 

Thus, using Eqs. (fT5|) and (|20p one can determine the 
integration constants A and B to be 



A 

where 



B 



e iSoL F-(L/2) + e- t5oL F+(L/2) 
APshxSoL 
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L 



L/2 



5{x)dx. 



(21) 



(22) 



-L/2 



It should be mentioned here that the solutions for 9 
and ip (Eqs. (ITU)) and ([T8j) . respectively) are not given 
explicitly but rather form a system of coupled equations. 
Fortunately, this system can be easily solved by a method 
of consecutive iterations. Indeed, starting with S(x) = 
one can solve Eq. (fTX?)) for 9 and next solve Eq. (fT5|) for 



C. Analytical approximation 

It is interesting that for the simplest, but widely 
used model of uniform bias current distribution ("f(x) = 
const), all the integrations can be performed analytically, 
leading to a compact and fully analytical expression for 
the I-V characteristic. 

Indeed, assuming the small-amplitude limit (tp <C 1) 
and ignoring self-pumping effects (a(x) — const), one 
can calculate the first-order approximation for Eq. (|11[) . 
Following the steps outlined above we find 9{x) — hx and 
consequently 



Ahx 



ip(x) 



, hsm[(h + S)L/2] iSx 

P(h 2 - S 2 ) P(h 2 - 5 2 )SsmSL 
hsm[(h-S)L/2] iSx 
P{h 2 -S 2 )5 sin SL 1 ' 



and 



70 = a r2+Im 



1 



h 2 (cos 5L — cos hL) 



2P{h 2 -S 2 ) LP(h 2 -~6 2 ) 2 Ssm5L 

(24) 

where P and <5 have been defined in Eq. (fTTj) . 

When the surface losses are neglected (/3 = 0) the 
above simple expression is reduced to the solution de- 
rived earlier—. Moreover, it can be shown that for (3 — 
the expression (|24|) is also equivalent to apparently dif- 
ferent solutions derived independently in Refs. Il5j|2l| in a 
form of infinite series expansions. 



D. Large-amplitude corrections 

As follows from Eq. (|TT|) . the evaluation of S(x) is cru- 
cial for the determination of I-V characteristic. So far 
we have assumed ^< 1 and used an approximation (1131) . 
However, for ip larger we should consider the exact rela- 
tion 



S(x) = (sin 0o cos-0)t + (cos 4>o sinip)- 



(25) 
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Since ip is an oscillatory function, one can use well- 
known relations 22 involving Bessel function. As shown 
in Ref . [HI, the leading terms for cos tp and sin tp are given 
by 



cosV' ^ M$\), 



smw ~ a. — —w, 



(26) 



where Jo, J\ denote the Bessel functions of order and 
1, respectively. 

Using Eq. (|2"6"|) we can find the time-independent con- 
tribution S(x) to bei& 



S( X ) = ^»lm 



ipe 



(27) 



Accordingly, the right-hand sides of Eqs. (fl4l and (fTo]) 
should be replaced by sin^ cos^ and Jo(\i>\)e 10 , respec- 
tively. 



E. Self-pumping effects 



According to Refs. Il6lll7l . the total current density due 
to the quasi-particle tunneling is given by 

oo 

Toff = J! J n (ev ac /huj) 7dc(«dc + nhu/e), (28) 



where J„ is the Bessel functions of order n, 7d c denotes 
the unpumped I-V dependence, and the total voltage 
applied to the junction can be separated into a constant 
(ud, c ) part and and an oscillatory part of amplitude v ac 
and frequency u). 

Using the Josephson relation^ t>dc = Tvuj/2e and com- 
ing back to our dimensionless notation we find v^ c = Q, 
w ac = \-tp\Q and 



7eff 



OO 

^J 2 (|^|/2) 7dc [r!(l + 2n)]. (29) 



For small arguments z = \ip\/2 <C 1 the infinite ex- 
pansion (|29[) can be truncated to include only quadratic 
terms in z: 



7cff = 7dc(«) + ^ [7dc(-fl) - 2 7dc (r!) + 7dc (3fi)] . (30) 

As an unpumped I-V characteristic we can take the 
nonlinear resistive model 2 ^ 



7dc(^) = a a Vl \ b 



(n/n g y 



W/n g )» + 1] 



(31) 



where n 3> 1, Q g denotes the normalized gap voltage and 
b = Rj/R n is the ratio of normal-state resistances below 



For n —¥ oo the highly nonlinear continuous depen- 
dence (|3ip tends to a simple discontinuous linear form 
which will be used in further calculations: 



7dc 



aQ$l(b 



1) 



for 



f2 < tig 

ft > fL 



(32) 



and above Q g . 



It is clear that 7 e ff is x-dependent due to the self- 
pumping via 4>{x), thus we can define an effective damp- 
ing factor a Q g(x) = ^ c q{x)/VL and compute it self- 
consistently by starting with a e g — olq, calculating ip, 
substituting into Eq. (|29j). evaluating a new approxima- 
tion for a e g(x), and so on. 

IV. RESULTS AND DISCUSSION 

In this section we compare analytical results derived 
above with numerical simulations obtained by the finite- 
difference implicit scheme^. The I-V characteristic is 
given by Eq. (fTTj) , while S(x) follows from the self- 
consistent solutions of Eqs. (flUl) and (|T5|) . To illustrate 
the influence of surface losses we consider first the sim- 
plest case 7 (x) = const where the I-V dependence can 
be expressed in a closed form (f2"4"|) . 

Fig. I shows the central part of the I-V characteristic 
for a junction of moderate length L — 5. The remain- 
ing parameters (an = 0.2, h = 3) are similar to those 
assumed in Ref. |l5|. Open circles denote the results of 
numerical simulations for a discrete set of 70 points. The 
dashed line denotes the analytical solution ([24|) while the 
solid line follows from the self-consistent solutions pH|) 

and (TTH1) obtained by consecutive iterations for 8 and ip. 
Fig. 1(a) shows the case f3 = (no surface losses) and in 
Fig. 1(b) we assume f3 = 0.01. Similar results for h = 5 
are presented in Fig. 2(a),(b). 

One can see that the I-V dependence departs from the 
Ohmic line 70 = aoQ in the region of ~ h, forming the 
main FF step modulated by a series of Fiske steps with 
voltage spacing Af2 ~ n/L. Analytical approximations 
are continuous and consist of a series of resonances while 
numerical results show typical hysteretic behavior and 
we observe only the segments of positive slope. 

It is clear that the presence of even very small surface 
losses changes significantly the I-V dependence. Gen- 
erally, the steps (resonances) become smaller and this 
effect is more pronounced for larger values of the exter- 
nal field h. Such a result can be easily explained if we 
recall that the main FF step and accompanying Fiske 
steps are visible in the region £1 ~ h and the surface loss 
factor enters the formalism via P = 1 + iCl/3. In other 
words, for higher external magnetic field the influence 
of surface losses is stronger and the I-V characteristic 
becomes more smooth. The influence of surface losses 
is clearly visible in experimental result a 16 ' 17 , where the 
Fiske steps gradually disappear as the external magnetic 
field is getting stronger. 

Comparing analytical and numerical results shown in 
Figs. 1 and 2 one can see that the fully analytical solu- 
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FIG. 1: (Color online) Current-voltage characteristic calcu- 
lated analytically (solid line) and numerically (open circles) 
for y(x) = const, L — 5, h — 3, ao = 0.2: (a) /3 = 0, (b) 
ft — 0.01. The dashed line shows the first-order approxima- 
tion pg|). 



tion (dashed line) is fairly accurate. On the other hand, 
the self-consistent solution (solid line) reproduces very 
accurately all the details of the numerical solution. 

As the next example, let us consider a more realistic 
case of a very long junction with small damping. Follow- 
ing Ref. [13 we choose L = 40, a = 0.033, (3 = 0.035, 
and a slightly asymmetric current profile 7(2) depicted 
in Fig. 3. The bias electrode (x\ < x < x%) is shorter 
than the total junction length. Consequently, the current 
distribution is assumed parabolic for x\ < x < X2, and 
exponential in the unbiased tails (x < x±, x > x%). 

Contrary to the previous example, now we cannot use 
the analytical approximation (|24p . First, the current dis- 
tribution is x-dependent, what means that 0(x) departs 
from a simple linear dependence 0(x) — hx, and rele- 
vant integrals cannot be calculated analytically. Second, 
it appears (even for a uniform current profile) that the 
first-order approximation is not accurate enough for very 
long junctions, and consecutive iterations are necessary 
to obtain a self-consistent solution. 
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FIG. 2: (Color online) the same as in Fig. 1 but for h = 5. 



Fig. 4 shows the I-V characteristics calculated for junc- 
tion parameters specified above and for two values of 
the external magnetic field h = 2.5 and h = 3.5. As 
before, open circles correspond to numerical simulations 
while the solid line represents a self-consistent solution. 
One can see that the agreement between numerical and 
analytical results is rather poor for h = 2.5. To ex- 
plain this discrepancy we should recall the main assump- 
tion (see Sec. II) that Eq. §5§ can be separated into a 
time-independent part and a part sinusoidal in fit while 
neglecting higher-order harmonics. The Fourier analy- 
sis shows, however, that the anharmonic contribution 
to tjj(t) is rather large for h — 2.5. E.g. for 70 = 0.3 
we find the content of the second harmonic to be about 
35%. Physically, it means that the fluxon train is not 
sufficiently dense for h = 2.5 and its time-dependence 
although periodic is not strictly sinusoidal. 

Contrary to the case h = 2.5, for h = 3.5 we observe ex- 
cellent agreement between numerical and analytical data. 
Now the I-V characteristic is smooth and the Fiske steps 
disappeared completely as a combined result of surface 
losses and self-pumping effects. The Fourier analysis of 
ip(t) shows now that the anharmonicity diminishes rather 
quickly with increasing h and for h = 3.5 we observe only 
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FIG. 3: Normalized bias current distribution similar to that 
assumed in Ref. 1 171 L = 40, x\ — —9, X2 = 5.5. 
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FIG. 4: (Color online) Current-voltage characteristic calcu- 
lated analytically (solid line) and numerically (open circles) 
for two values of the external magnetic field: h = 2.5 and 
h — 3.5. The current profile is shown in Fig. 3 and the remain- 
ing junction parameters are: L = 40, ao = 0.033, /3 = 0.035. 



about 10% of the second harmonic. 

Another interesting feature which is visible for h = 3.5 
is an additional small step at Oj, = 1.9. Such a step ap- 
pears in experimental result a 17 i 26 and can be attributed 
to the self-pumping effect described in Subsection III.E. 
As shown in Ref. |2g, the position of the current step fol- 
lows from a simple relation Qb = £l g /3, where Q g denotes 
the gap voltage. Assuming typical junction parameters 26 
we find a normalized dimensionless gap voltage Q g = 5.7, 
hence = 1.9. 

It is interesting that the value of Of, can also be 
deduced directly form Eq. (|30l) . Indeed, noting that 
7d c (— ^) = — 7dc(^) one can see the quadratic term in 
Eq. pop to vanish, provided 7dc(^) is a linear function. 
However, for a nonlinearity starting at £1 = £l g we find 
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FIG. 5: (Color online) Spatial distribution of a e a(x) calcu- 
lated analytically (solid line) and numerically (dotted line) 
for ft = 3.5 and two values of the averaged current density 
70 = 0.1 and 70 = 0.3. The remaining junction parameters 
are the same as in Fig. 4. 



7dc(3f2) 7^ 37dc(^), thus a nonzero contribution to 7 c g 
appears at 30 = fi s , giving rise to a step at f2 b = fi g /3. 

As shown in Subsection III.E, an effective loss factor 
a e g(x) is spatially modulated due to the self-pumping, 
and its average value becomes significantly larger than 
ao- Fig. 5 shows the spatial distribution of cv c ff plotted 
for 70 = 0.1 and 70 = 0.3. The junction parameters are 
assumed as before, i.e. L = 40, h = 3.5, a = 0.033, 
P = 0.035. Now the dotted lines correspond to numeri- 
cal simulations and the solid lines denote self-consistent 
(analytical) solutions. One can see that the agreement is 
excellent, the self-consistent solution following closely all 
the details of the numerical simulations, taken here as a 
reference. 

It is clear that a e g generally grows with increasing 
value of 70. E.g. for 70 = 0.3 we obtain an averaged 
value of a e f{(x) about three times larger than the "un- 
pumped" value ao- Such an effect together with surface 
losses makes the I-V curve smooth, damping effectively 
the Fiske steps. 



V. CONCLUSIONS 

In this paper a semi-analytical approach has been sug- 
gested, making possible to solve a modified sG equa- 
tion ^ with both surface losses and self-pumping effects 
taken into account. The solution, as given by Eqs. Q, 
(fTU|) . and (fl"8|) . consists of a rotating background and a 
unidirectional fluxon train accompanied by two plasma 
waves traveling in opposite directions with the velocity 
close to the critical value £l/So ~ ±1. Having obtained 
an analytical solution to Eq. ([2]) one can determine some 
macroscopic, directly measurable quantities, such as the 
constant bias current and constant voltage, making pos- 



7 



sible to calculate the I-V characteristic (fTTj) of the junc- 
tion. 

For the uniform current density distribution (7(2;) = 
const) it has been shown that the relevant expressions for 

6(x), ip(x), and consequently the I-V dependence can be 
obtained in a closed fully analytical form (j2~4")l . As fol- 
lows from Figs. 1 and 2, such an approximation is fairly 
accurate for moderate junction parameters. In general, 
however, practical FF oscillators are based on very long 
junctions with strongly nonuniform current profile. In 
such a case an analytical approximation (|24p appears in- 
sufficient and one should look for a self-consistent solu- 
tion of Eqs. (fT0|) and (fTS)) which can be obtained by an 
iterative procedure. 

As shown in Figs. 1,2,4 and 5, the self-consistent so- 
lutions and numerical simulations are generally in very 
good agreement. The only exception, where one can ob- 
serve qualitative rather than quantitative agreement is 
the I-V curve shown in Fig. 4 for h — 2.5. As explained in 
the previous Section, the fluxon train is not dense enough 
for h = 2.5 and consequently the time-dependent part of 



the solution is periodic but not strictly sinusoidal, vio- 
lating the main assumption of the present approach. 

Finally, it should be mentioned that the method pre- 
sented here can be easily extended to include general, 
more realistic boundary conditions discussed in Refs. 
[HUTt]!!?]. In the present paper, however, we intention- 
ally assumed standard open-circuit boundary conditions 
([3]), where the plasma waves could interfere after com- 
plete reflection at the boundaries, giving rise to clearly 
visible Fiske steps. This way, we were able to separate the 
influence of surface losses and self-pumping from that of 
boundary conditions, which (e.g. for a resistive load) can 
also affect the I-V characteristic and make the physical 
mechanism discussed here less clear. 
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